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STEADY- STATE  POINT- SOURCE  EXCITATION  OF  A  LAMINATED  COMPOSITE 

Abstract 

Steady-state  periodic  excitation  at  a  point  of  an  extended  periodically 
laminated  elastic  composite  Is  considered  In  anti-plane  strain.  The  curves 
of  constant  phase  are  determined  In  the  geometric  optics  approximation.  The 
associated  distribution  of  group  velocity  is  also  calculated. 
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Introduction 

The  propagation  of  horizontally  polarized  harmonic  shear  waves  through 
a  periodically  layered  elastic  body  of  unbounded  extent  Is  examined  here  In 
the  geometric  optics  approximation*  The  geometric  optics  approach  Is  valid 
for  high  frequencies  where  the  relevant  wavelength  Is  small  In  comparison 
with  some  characteristic  dimension  of  the  problem*  Under  this  approximation 
any  quantity  which  describes  the  wave  field  can  be  expressed  In  the  form  of 
an  amplitude,  which  varies  slowly  In  space  and  time,  multiplied  by  a  phase 
factor  of  the  form  e  ,  where  0  Is  the  frequency  and  S  Is  the  phase* 
Following  the  notation  of  reference  [1],  the  phase  S  can  be  written  In  the 
form  S  -  t-o(x)*  The  surfaces  S  ■  constant,  which  are  surfaces  of  constant 
phase  will  be  examined  for  the  steady  state  case  of  a  constant  frequency 
point  source* 

Infinite  Medium 

The  propagation  of  horizontally  polarized  harmonic  shear  waves  through 
a  periodically  layered  body  of  Infinite  extent  has  been  analyzed  In  some 
detail  in  reference  [2]  and  is  summarized  here  as  far  as  needed  In  the 
following* 

r 

The  system  to  be  considered  consists  of  an  infinite  sequence  of  two 
alternating  layers,  each  of  which  is  homogeneous,  isotropic  and  perfectly 
bonded  to  the  adjoining  layers*  A  unit  cell  is  defined  as  the  union  of  any 
two  adjacent  layers*  The  two  lamellae  of  a  typical  unit  cell  have  elastic 
constants  X,  p;  X  ,  v  ,  thicknesses  2h;  2h',  and  densities  p;  p',  see 
Fig*  la* 


By  utilizing  the  equations  of  elasticity  and  Imposing  continuity  con¬ 
ditions,  which  specify  that  the  displacements  and  stresses  are  continuous  at 
each  layer  Interface,  we  can  write  the  equation  of  motion  for  the  laminated 
body.  For  the  Infinite  medium  this  equation  takes  on  the  form  of  a  partial 
differential  equation  with  coefficients  which  are  periodic  In  the  direction 
perpendicular  to  the  layering  with  period  d  »  2h  +  2h''.  Application  of  the 
one-dlmenslonal  theory  of  Floquet  for  differential  equations  with  periodic 
coefficients  leads  to  the  result  that  the  details  of  the  motion  of  the 
entire  laminated  body  can  be  obtained  by  considering  only  one  unit  cell. 
Therefore  the  boundary  conditions  for  one  unit  cell  using  Floquet 's  theory 
represent  a  set  of  four  homogeneous  equations  which  are  sufficient  to  de¬ 
scribe  the  entire  laminated  body. 

For  a  nontrivial  solution  the  determinant  of  the  matrix  of  coefficients 

of  these  four  homogeneous  equations  must  vanish,  yielding  the  dispersion 

equation  H'  (i2,  rt*  ~  0  In  expanded  form 

2 

^yaa'  cos  +  e)  +  (yo  “a")  cos  x(a  ~  ea") 

“(ya  +  a')^  cos  irCa  ea')  ■  0  (1) 

V2  2  I  2  2  2 

Q  ~C.a'*Vo  fi  -  C.  where  fi,  5,  and  t\  are,  respectively, 
the  nondimensional  frequencies  and  wave  numbers  defined  by 


a 


(2) 


where  u,  k^,  and  k^  are,  respectively,  the  circular  frequency  In  radians 
per  unit  of  time,  the  wave  number  in  the  direction  parallel  to  the  layering. 


and  the  wave  number  In  the  direction  perpendicular  to  the  layering.  In 
addition,  the  geometric  and  material  parameters  are  defined  in  the  following 
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manner:  y  Is  the  ratio  of  shear  moduli  (p/y'),  e  Is  the  ratio  of  layer 
thicknesses  (2h'’/2h),  and  a  is  the  ratio  of  shear  wave  speeds 

o  *  /pp'/p'p 

Since  the  dispersion  equation  related  a  nondlmenslonal  frequency  to 
two  nondlmenslonal  wave  numbers,  the  roots  of  this  equation  define  a  surface 
in  three-dimensional  frequency-wave  number  space.  A  qualitative  sketch  of 
this  surface  is  shown  in  Fig.  lb.  It  is  seen  that  the  surface  is  discon¬ 
tinuous  at  n  ~  n/(l  +  e)>  n  «  1,  2,  3,...  These  planes  of  discontinuity 
divide  the  surface  into  bands  known  as  Brlllouln  zones.  At  the  end  of  each 
Brillouin  zone  (with  the  exception  of  n  *  0)  the  dispersion  equation  admits 
complex  as  well  as  real  roots.  These  portions  of  the  dispersion  spectrum 
for  which  the  wave  numbers  are  complex  are  called  stopping  bands  and  corre¬ 
spond  to  regions  of  the  spectrum  where  harmonic  waves  are  not  propagated 
but  attenuated  exponentially.  Curves  of  constant  frequency  on  the  surface 
over  the  first  two  Brlllouln  zones  are  shown  in  Fig.  2.  These  curves  play 
an  Important  part  in  the  construction  of  the  surfaces  of  constant  phase. 

Geometric  Optics 

Introducing  the  geometric  optics  approximation  we  assume  that  any 
quantity  which  describes  the  wave  field  can  be  expressed  by  a  formula  of 
the  type  Ae  where  A  is  the  amplitude,  which  varies  slowly  in  space 
and  time,  (2  is  the  frequency,  and  S  is  the  phase.  Following  the  notation 
and  methodology  outlined  in  reference  [1],  we  write  S  in  the  form 
S  ■  t  -  a(*)*  The  wavefronts  defined  by  S  (x,  t)  “  0  can  be  found  by 
examining  the  surfaces  a(x)  -  constant  for  fixed  time.  These  "snapshots” 
give  the  successive  positions  of  the  wavefront  in  physical  space.  To  find 


the  shape  of  these  wavefronts  In  physical  space  we  begin  by  rewriting 
the  dispersion  equation  for  the  steady  state  case  of  a  single  frequency, 
say  Q  ■  Qq.  Doing  this  the  dispersion  relation  takes  the  fora 

H  (Pl.  P2)  “0  pj  -  (3) 

where  the  p^  can  be  Interpreted  as  the  components  of  the  slowness  vector 
defined  by  dispersion  relation  In  this  fora  Is  called 

the  elkonal  equation  or  equation  of  geometrical  optics.  The  standard 
method  of  solving  equation  (3)  Is  by  means  of  the  characteristic  or  ray 
equations,  which  proceeds  as  follows,  ref.  1: 

Consider  an  equation  In  n  Independent  variables  (x^,  ... 

and  a  function  g  (x^,  x^,  x^,  ...  x^)  which  satisfies  a  differential 
equation  H  (p,  0,  x)  where  p  and  x  are  vectors  of  dimension  n  and 

1-1.  2.  3....  (4) 

We  wish  to  find  whether  there  are  any  curves  In  x-space  for  which  the 
differential  equation  H  (p,  a»  x)  reduces  to  a  set  of  first  order  ordinary 

differential  equations. 

We  begin  by  noting  that  any  curve  In  x-space  may  be  written  In  the 
parametric  form  x  ■  x  (x).  The  total  derivative  of  o  along  one  of  these 

curves  Is  given  by 
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where  the  summation  convention  Is  employed.  In  addition  consider  the  total 
derivative  of  along  the  same  curve 

2 

dPj  d  3o  3  o  dx 

■  m  —  -  f  '\  m  —  -  -  J 

dx  dx  '■axj-'  ax^axj  dx 

Finally  let  us  take  the  x^  derivative  of  H  (p,  o,  x)  ■  0  which  gives 


JLJ—  a^L+iH 

3*j3*j  3Pj  3o  3*i  3*^ 

By  choosing  as  our  definition  for  the  special  curves  In  x-space  the 
following 


(7) 


dXj  3H 
^  "  3Pj 

we  may  combine  equations  (6)  and  (7)  to  give 


(8) 


dp^  3H  3H 

Tx - •  111 

The  use  of  equation  (7)  In  (4)  then  yields 


(9) 


ill 

dx 


(10) 


The  set  of  equations  (8),  (9)  and  (10)  now  represents  a  complete  set  of 
(2  n  +  1)  ordinary  differential  equations  for  determining  a  characteristic 
curve  xj  (x)  and  the  values  of  a  and  p^  along  It. 


Specializing  these  results  to  our  case  gives 


dx.  dp  da 

■dr=(V^i'  ^“0*  ^ 


(P-^  H) 


(11) 


«diere  A  is  a  parameter  along  the  characteristic  curve  (ray).  The  second  of 
equations  (11)  Indicates  that  pj  are  constant  along  rays,  hence  the  ray 
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direction  VpH  is  constant  and  therefore  the  rays  are  straight  lines.  Since 

the  pji  are  constant  we  can  integrate  equation  (11)  directly.  For  the  case 
of  a  point  source  of  constant  frequency  located  at  the  origin,  we  can  readily 
obtain, 

X.  -  i.s  o  ■  P«J18  (12) 


where  £  is  a  unit  vector  defined  by 


I 


(13) 


and  8  is  the  distance  from  the  source  measured  along  the  ray.  From  these 
equations  we  find  that  the  wavefront  S  =  0  is  Ipcated  at 


s 


(14) 


and  the  coordinates  of  a  point  on  the  wavefront  are 


(15) 


By  varying  pj  over  all  values  satisfying  H  (p^,  P2)  “  0  for  fi  ■  we  can  find 

the  entire  wavefront.  It  should  be  noted  that  these  wavefronts  can  be  found  by 
examining  curves  of  constant  frequency  in  p-space  (Fig.  2).  If  for  example  we 
choose  a  particular  value  for  say  0.9,  (see  Fig.  3),  then  i  is  merely  the 


unit  normal  to  the  curve  for  any  choice  of  p^  satisfying  H  (p^,  P2)  *  0  fot 
n  ■  0.9.  By  varying  pj^  along  the  curve  we  can  trace  out  the  entire  wavefront 
(Fig.  4a). 


With  the  aid  of  a  compass  and  two  triangles  a  very  simple  geometric  con¬ 
struction  of  these  wavefronts  Is  possible.  We  begin  by  drawing  a  coordinate 


system  with  the  sane  orientation  as  the  coordinate  system  In  Fig.  2  From 
Fig.  2  choose  a  value  for  and  select  a  point  along  the  curve  n  >  constant. 
Now  draw  the  normal  to  the  curve  at  this  point.  With  the  aid  of  the 
triangles  draw  a  line  passing  through  the  origin  of  our  new  coordinate 
system  which  Is  parallel  to  this  normal.  Now  working  with  Fig.  2  project 
the  wave  number  vector  onto  the  unit  normal  vector.  The  length  of  this 
vector  Is  proportional  to  p.i  (since  g,  is  a  unit  vector).  All  that 
remains  Is  to  take  the  Inverse  of  this  length.  This  is  accomplished  by 
drawing  a  reference  hyperbola  and  using  the  compass  to  measure  the  Inverse. 
This  length  Is  then  measured  from  the  origin  of  our  new  coordinate  system 
along  the  paralled  to  g.  This  Is  the  first  point  of  the  wavefront.  The 
process  Is  repeated  for  different  points  on  the  same  curve  n  *  constant 
until  a  sufficient  number  of  points  exist  In  order  to  pass  a,  smooth  curve 
through  them. 

Because  the  elkonal  equation  Is  unchanged  by  switching  from  p^  to 
-p£  the  wavefronts  will  be  symmetric  in  the  X2  axis.  The  complete 
picture  Is  shown  In  Fig.  5.  The  wavefronts  for  different  values  of  are 
shown  In  Fig.  4b,  4c,  4d. 

These  figures  of  course  give  no  Information  concerning  energy  dis¬ 
tribution.  It  may  be  of  Interest , though,  to  examine  the  magnitude  and 
direction  of  velocity  of  energy  propagation.  To  accomplish  this  a 
numerical  finite  difference  calculation  was  carried  out  to  show  the 
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group  velocity  vector  at  points  on  the  zero  phase  curves.  Here  the 
group  velocity  Is  defined  as; 


where  H'ls  the  dispersion  relation  written  as  a  function  of  ii,  n 
which  were  previously  defined  and  e^,  e^  which  are  unit  vectors  In  the 
direction  perpendicular  and  parallel  to  the  layering, respectively. 

The  Figs, 6a,  6b,  6c,  6d  duplicate  the  phase  curves  previously  shown .but 
now  Include  group  velocity  Information  as  well.  The  magnitude  of  the 
group  velocity  vector  normalized  by  the  phase  velocity  /p/p  '  Is  drawn 
radially  outward  from  each  zero  phase  curve  where  the  vectors  begin,  not 
at  the  origin  but  at  the  appropriate  point  on  each  phase  curve.  It 
should  be  noted  that  the  zero  phase  curve  passing  through  the  conical 
point  (Fig.  6c)  has  a  polntwlse  discontinuity  In  the  value  of  the  group 
velocity  at  the  conical  point. 
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Conclusion 

Through  the  utilization  of  the  geometric  optics  approximation  surfaces 
of  constant  phase  were  examined  for  the  case  of  a  point  source  Imbedded  In 
a  periodically  layered  elastic  body.  Knowledge  of  the  shape  of  these 
surfaces  Is  the  first  step  toward  understanding  how  more  complicated  signals 
propagate.  For  Instance,  we  already  know  from  examining  these  zero  phase 
surfaces  that  an  observer  stationed  somewhere  within  the  medium  would  detect 
one,  two,  or  three  distinct  signals  depending  upon  ones  orientation  with 
respect  to  the  source.  A  similar  phenomenon  occurs  In  the  case  described 
In  reference  [1]  for  waves  In  a  MHD  plasma. 
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Fig.  3  Curve  of  H  ■  0.9  showing  vectors  £  and  £. 

Fig.  4a,  b.  Curves  of  constant  phase  for  different  values  of  Q. 
c,  d 

Fig.  5  Complete  wavefront  for  0  ~  1.0. 

Fig.  6a,  b.  Distribution  of  group  velocity  on  curves  of  constant  phase, 
c,  d 
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Figure  6a 


